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Abstract 

To be able to construct viscoelastic material models from fractional-order differin- 
tegral equations that are applicable for 3D finite-strain analyses requires definitions for 
fractional derivatives and integrals of symmetric tensor fields, like stress and strain. 
Here we define these fields in the body manifold. We then map them into spatial fields, 
expressed in terms of an Eulerian or Lagrangian reference frame where most analysts 
prefer to solve boundary value problems. 


1 Definitions 


Liouville and Riemann 1 defined fractional- order integration as an analytic continuation of 
Cauchy’s n-fold integral by writing [4, f5, Eqn. A] 

ry(x):= W)L W^W^ yix ' )dx '’ a ’* eI4 ’ (1) 

where J a is the Riemann-Liouville integral operator of order a. 

From this single definition for fractional integration, one can construct several definitions 
for fractional differentiation (cf. e.gi, [10, 13]). The special operator D% that we choose to 
use, which requires the dependent variable y to be continuous and [a] -times differentiable 
in the independent variable x, is defined by 

D; y(x) := jW-»DWy(j), (2) 


such that 


D+J a y(x) — y(x) and J a D*y(x) = y(x) 


H k 

Z_y k \ y 0+ 


k=0 


a 


6 M+, 


(3) 
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with yjp := D*y(0 + ) wherein D n , n G N, is the classical differential operator. It is accepted 
practice to call D% the Caputo 2 differential operator of order a, since Caputo [1] was the 
amoung the first to use this operator in applications and to study some of its properties. 


2 Continua 


Body-tensor fields, space-tensor fields, Cartesian space-tensor fields, and the mappings be- 
tween them have all been carefully documented by Lodge in [5, 6, 7]. 

A continuum consists of an infinite set of point particles {*}3}, also referred to as mass 
elements, filling a region B in 3-space (i.e., ® C R 3 ). We call this set the body B. In any 
admissible body-coordinate system B , defined over ®, each particle ^3 in B is assigned a 
unique set of body coordinates, { — (^ 1 ,^ 2 ,^ 3 ), \ l 6 R, that are independent of time (i.e., 
B: *}3 — )• <£, cf. Lodge [6]); they convect with the body. 

This same continuum ® can also be represented by an infinite set of point places {3^o } 
occupying a connected region in space S at some arbitrary time to denoting its reference 
state. Each place 3£o relates to a unique particle *}3 in B and is given a label of X, which 
corresponds to the spatial position of 3£o (and therefore of ^3) in this reference configuration. 
Given an admissible, rectangular-Cartesian, coordinate system C defined over S, each place 
3£o in S is thereby assigned a unique set of spatial coordinates X = (Xi,X 2 ,X 3 ), X* 6 R, 
that are fixed in space (i.e., C: X o — » X). 

Later, at current time t ( t > to), continuum B coincides with another infinite set of 
point places {3f} that now occupies a different region in space S. Each place X relates to a 
unique particle in B and is given a label of x with coordinates x — (xi, X2, X3) in C, which 
corresponds to the spatial position of X (and therefore of ^3) in this current configuration. 
Because C: X 0 — > X and C : X — >• x, the location of a particle in space varies over time. 

Particle ^3 moves through space § with a velocity v(t) of 


v 


dx 

Ik' 


( 4 ) 


whose components — dxi/dt, i = 1, 2, 3, are quantified in the Cartesian coordinate system 
C. 

The fundamental hypothesis of Cartesian continuum mechanics is that the motion of 
each particle in the body is assumed to be sufficiently smooth in the sense that mappings 


Sx = F-5X and Sv = DF ■ SX = DF ■ F~ l ■ dx = L ■ Sx (5) 


exist, where F(to,t) := dx/dX defines the deformation- gradient tensor and L{t) := dvjdx 
the velocity-gradient tensor, neither of which is symmetric. The deformation gradient F is 
positive definite because, from the conservation of mass, 

0 < — det F < 00 , (6) 

2 Actually, Liouville introduced the operator in his historic first paper on the topic [4, ^[6, Eqn. B]. 
Still, nothing in Liouville’s works suggests that he ever saw any difference between D+ = j^ a ^~ a D^ a ^ and 
D a = D a being his accepted definition [4, first formula on pg. 10] — the Riemann-Liouville 

differential operator of order a. Liouville freely interchanged the order of integration and differentiation, 
because the class of problems that he was interested in happened to be a class where such an interchange 
is legal, and he made only a few terse remarks about the general requirements on the class of functions for 
which his fractional calculus works [8]. Rabotnov [11, pg. 129] also introduced this differential operator into 
the Russian viscoelastic literature a year before Caputo’s paper was published. 
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and consequently F 1 (^o 5 ^) — dX/dx exists. In contrast, L is not positive definite, in 
general, and as such L~ l need not exist. Here g = g(t) denotes mass density with po — Q{to)- 

3 Field Transfer: Body To Cartesian Space 

The operation of field transfer makes it very plain as to whether a particular spatial field 
is Eulerian or Lagrangian; it is a consequence of the time when field transfer takes place. 
Eulerian fields result from a transfer of field at current time t given by the mapping: body 

field t==£- space field ; whereas, Lagrangian fields result from a transfer of field at reference 

to 

time to (which is often taken to be zero) given by: body field t=> space field. Details of the 

mathematics underlying the mappings t=> and can be found in the texts of Lodge [5, 7]. 
An important property of the field transfer operator is that the resulting spatial fields are 
objective (i.e., frame invariant). 

Let us consider (i) an arbitrary, symmetric, covariant, body tensor and (ii) an 

arbitrary, symmetric, contravariant, body tensor 77 (^ 3 ; t), whose mappings into Cartesian 
space are known; specifically, let the covariant field map as 




<0 


M(X;t), 

N(X 0 ;t 0 ,t), 


such that N = F t - M • F, 


(7a) 


and the contravariant field map as 

t 

to 


v PP; t) 


C?(3£;«), 




such that H — F~ l - G ■ F ~ T , 


(7b) 


where M and G are some arbitrary (but known), symmetric, Eulerian, tensor fields, with N 
and H designating their respective, symmetric, Lagrangian counterparts. In these transfor- 
mations of field, the deformation gradient F serves as a Jacobian of transformation between 
the two Cartesian frames that pulls the known Eulerian field backwards, out of the Eulerian 
frame and into the Lagrangian frame. 

3.1 Rates of Tensor Fields 

Lodge [5, pp. 321-327] has shown that the time rate-of-change D (= d/dt ) of covariant 
tensor pi maps into Cartesian space as 


Dpi ’ suc h DN = F T ■ M ■ F, 

° DN, 


(8a) 


while the time rate-of-change of contravariant tensor rj maps as 

t 


Dr] 


to 


G ’ such that DH = F~ l -G-F~ T , 
DH , 


(8b) 


wherein 


M := DM + v • VM + L T ■ M + M • L and G DG + v VG - L ■ G - G L T ( 9) 
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denote the lower- and upper- convected derivatives, respectively, of Oldroyd [9], which reduce 
to Lie derivatives taken with respect to velocity v whenever DM — 0 or DG = 0. The 
vector operator V represents the spatial gradient d/dx. 

Even though Cartesian tensor fields cannot distinguish between covariance and con- 
travariance, rates-of-change of Cartesian fields in the Eulerian frame do have an intrinsic 
dependence upon this property. 


3.2 Fractional Rates of Tensor Fields 

The Caputo derivative D+ (defined in Eq. 2) of covariant tensor /x maps into Cartesian 
space as 

{ t=L. D aV M 

tn * ’ such that D?N = F t - D; v M • F, (10a) 

^ D?N, 

while the Caputo derivative of contravariant tensor rj maps as 


d; n 


to 


d; a g, 

d;h. 


such that 


D°H = F 


D? G ■ F 


(10b) 


wherein 

D^M := 


T(1 


1 f 1 1 

- <x) Jto (t - 


F~ T (t', t) ■ M(t') • F~ l {t' , t) dt', 0 < a < 1, (11a) 


t>y 


and 


0<a<1 ’ (llb) 

are objective rates of order a, with operator D* v being affiliated with covariant- like fields 
and D+ A with contravariant-like fields. Unlike D+N and D+H, which are actual Caputo 
derivatives, derivatives D “ v M and D+ A G are not true Caputo derivatives, in a strict sense 
of the definition, which is why they are given different notations. 

A rigorous derivation of these results is given in Ref. [3]. 


3.3 Fractional Integrals of Tensor Fields 

The Riemann-Liouville integral J a (defined in Eq. 1) of covariant tensor /x maps into Carte- 
sian space as 


J aV M 

J a li u ’ such that J a N = F T - J aX7 M ■ F, 


J a N , 

while the Riemann-Liouville integral of contravariant tensor rj maps as 


(12a) 


r v 


to 


J aA G, 

J a H, 


such that J a H - F~ x ■ J aA G ■ F~ T 


(12b) 
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The Eulerian tensors J aV M and J aA G are objective integrals of order a defined by 


J aV M := -fi- [* ■ M(t') ■ F-'(t\t)d£ (13a) 

r(a) Jto {t - t'y 

and t 

jaA ° = W)L ¥^ F(t ' t} ■ G(0 • FT(t ' t)dt ’ (i3b) 

Unlike J a N and J a H , which are actual Riemann-Liouville integrals, integrals J qV M and 
J aA G are not true Riemann-Liouville integrals, in a strict sense of the definition, which is 
why they are given different notations. 

Again, a rigorous derivation of these results is given in Ref. [3]. 


4 Some Viscoelastic Models 


To illustrate the process of analytically continuing a known constitutive equation to one of 
fractional order, we consider three examples: the rubberlike liquid, the convected Maxwell 
fluid, and the convected Kelvin- Zener solid (cf. Lodge [5]). Expressed in terms of body 
tensor fields, the rubberlike liquid is 


7v(t) + p 7 1 (t) = fl 





the convected Maxwell fluid is 


r Dir(t) + 7r(t) + p 7 1 (t) = -r) Dj : (t), 


(14a) 


(14b) 


and the convected Kelvin-Zener solid is 

t Dn(t) + 7 T (t) + P7 _1 W = M 7 _ 1 (to) - V - D 7 _1 ( i ) 5 ( 14c ) 


where n is the symmetric, contravariant, stress tensor, 7 is the symmetric, 
postive-definite, covariant, metric tensor, M is the memory function (usually taken to be a 
decaying exponential), and r, 77 and /j, are material constants. All three of these models are 
assumed to be incompressible, with p being a Lagrange multiplier introduced to force this 
constraint: det 7 (t) = det 7 (to)- 

In much the same way that Caputo and Mainardi [2] analytically continued models in ID, 
we analytically continue the 3D models in Eq. (14) by writing the fractional-order rubberlike 
liquid as 

7z(t) + P7 _1 (t) = fi J a 7 _ 1 (t), (15a) 

the fractional-order convected Maxwell fluid as 


T a D?n(t) + 7 r(t) + p 7 ~\t) = W, (15b) 

and the fractional-order convected Kelvin-Zener solid as 

r Q D“ 7r(t) + 7T (t) + P7” 1 (0 = ^ 7 _ 1 (^o) - V a D*'y~ 1 (t). (15c) 

The memory function in Eq. (15a) is now considered to be an Abel (power law) kernel 
instead of a Boltzmann (exponential) kernel. The fact that the fractional derivatives on 
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both sides of the equation in Eqs. (15b & 15c) are of the same order (viz., a, 0 < a < 1) 
ensures that stress waves propagate with finite velocities in these models, which is important 
both physically and numerically. 

In an Eulerian transfer of field, Lodge [5, 7] has shown that 7r (t) T(t ), 7 _1 (t) i=^ /, 

t t A 

7 -1 (t 0 ) )=>- B(to,t), and I — —2 D(t), where T is the Cauchy stress tensor, I 

is the identity tensor, B (:= F-F T ) is the Finger deformation tensor, and D (:= + L T )) 

is the rate-of-deformation tensor. Using these results, along with Eqs. (10 & 12), one can 

apply the field transfer operator t===> to Eq. (15) and get the following fractional-order 
viscoelastic models in the Eulerian frame. The fractional-order rubberlike liquid becomes 

r + pI = ,r> I = JL- ) £ Bit', t ) < (16a) 

the fractional- order upper-convected Maxwell fluid becomes 

t *DZ a T +T+pI= 2r] a D; A D , (16b) 

and the fractional-order upper-convected Kelvin-Zener solid becomes 

r a D« A T + T + pI = fiB + 2 rj a D? A D, (16c) 

where now the constraint for incompressiblity is det F = 1 in accordance with Eq. (6). Their 
Lagrangian counterparts can be written straightaway, if needed. 
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